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We report a detailed investigation of GdOi-^F^FeAs (x = 0, 0.07 and 0.14) samples by means of 
high- field/high- frequency electron spin resonance (HF-ESR) together with measurements of thermo- 
dynamic and transport properties. The parent GdOFeAs compound exhibits Fe long-range magnetic 
order below 128 K, whereas both doped samples do not show such order and are superconducting 
with T c = 20 K (x = 0.07) and T c = 45 K (x = 0.14). The Gd 3+ HF-ESR reveals an appreciable 
exchange coupling between Gd and Fe moments, through which the static magnetic order is clearly 
seen in the parent compound. Owing to this coupling, HF-ESR can probe sensitively the evolution 
of the magnetism in the FeAs planes upon F doping. It is found that in both superconducting sam- 
ples, where the Fe long-range order is absent, there are short-range, static on the ESR time scale 
magnetic correlations between Fe spins. Their occurrence on a large doping scale may be indicative 
of the ground states' coexistence. 



I. INTRODUCTION 

Iron-pnictide superconductors 1 - with superconducting 
critical temperatures up to 55K2r— have attracted a 
huge interest due to striking similarities to supercon- 
ducting cuprates as well as due to their original prop- 
erties. Indeed, most families of these layered materials 
feature an antiferromagnetically (AFM) ordered parent 
compound, and the evolution of superconductivity con- 
comitantly with suppression of AFM order upon doping. 
However there are important differences which render the 
Fe-pnictidcs a separate new class of superconducting ma- 
terials. Most striking of them are semi-mctallicity and 
the spin density wave (SDW) character of the AFM or- 
der in the undoped pnictides contrasted with the Mott- 
insulating AFM state in the cuprates, as well as a multi- 
band versus single-band electronic structure in the Fe- 
pnictide and cuprate high-temperature superconductors, 
respectively. 

Beyond study of the superconducting ground state, 
and of the magnetic and associated structural transitions 
seen in the parent compound, much attention has been 
devoted to the issue of the ground states' coexistence. 
Discrepancies on this issue have been found between dif- 
ferent families^ - — , with the variation of the boundary 
of the two ground states and different length scales of co- 
existence, especially in the so-called 1111 family. In this 
family, which has the composition ROFcPn (R - rare 
earth, Pn - pnictide) , the superconductivity evolves with 
the substitution of fluorine for oxygen. Here, replace- 
ment of one rare earth element with another can cause a 
significant variation of properties. Whereas in La-based 
superconducting samples there is evidence against static 
magnetic order in the FeAs planes^, in the case of su- 
perconducting samples based on different magnetic rare 



earths (R = Sm, Nd, Ce) evidence of remanent static 
magnetism is founc&ii. The situation appears compli- 
cated due to the fact that the magnetism then tends to 
be of a short-range order or disordered, possibly even 
dynamic^, which calls for the use of local probe tech- 
niques. These two different pictures complicate the es- 
tablishment of the unified phase diagram for 1111 pnic- 
tides, necessary for the full understanding of these mate- 
rials. In addition, as was shown by NM R 13 i 14 and /xSRp£ 
studies, there is a magnetic coupling between 4/ (Ce, Pr 
and Sm) and 3d (Fe) moments. Such coupling of the rare 
earth to the FeAs plane might give an additional contri- 
bution to the difference in physical properties of different 
1111-type superconductors. 



In the present work we investigate the evolution of the 
magnetism upon fluorine doping in Gd-based 1111 com- 
pound by means of high field/high frequency electron 
spin resonance (HF-ESR) complemented with measure- 
ments of thermodynamic and transport properties. The 
ESR data reveal a significant exchange coupling of Gd- 
and Fe-moments in the parent GdOFeAs sample which 
enables the Gd 3+ HF-ESR to probe sensitively the for- 
mation of the static SDW magnetic order in the FeAs 
planes. Interestingly, it is found that the signatures of 
such an order are still observed in the ESR spectra af- 
ter doping. In particular, though long-range SDW order 
present at very low doping is suppressed at doping lev- 
els where superconductivity appears, our results imply 
static on the ESR time scale, likely short-range, magnetic 
correlations between Fe spins. This result suggests that 
GdOi-^F^FeAs compounds may feature coexistence of 
quasi-static magnetism and superconductivity on a large 
doping range. 
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II. EXPERIMENTAL 



A. Setups 



The magnetization has been studied by means of a 
commercial SQUID magnetometer (MPMS-XL5, Quan- 
tum Design). For the thermal expansion measurements 
a capacitance dilatometer was utilized, which allows a 
very accurate study of crystal length changes. We mea- 
sured the macroscopic length changes dL/L of polycrys- 
talline samples. The linear thermal expansion coeffi- 
cient a was calculated as the first temperature deriva- 
tive of dL/L, while the volume thermal expansion co- 
efficient is given by /3 vo i = 3a for our polycrystallinc 
samples. The specific heat was studied in a Quantum 
Design PPMS calorimeter by means of a relaxation tech- 
nique. In the electrical transport experiments the sam- 
ples were investigated by four-probe p measurements us- 
ing an alternating DC-current. The ESR measurements 
at a frequency of v = 9.6 GHz were carried out in a 
standard Brukcr EMX system. The HF-ESR experi- 
ments were performed with a home-made spectrometer— 
at frequencies v = 83 — 348 GHz and magnetic fields 
B = — 15 T. All ESR measurements were made in a 
temperature range of 5 — 300 K. 



B. Sample preparation 

The polycrystalline samples GdOi-^F^FeAs (x = 
0,0.15,0.17, nominal content) were prepared by two dif- 
ferent routes. Route 1, which is similar to that described 
in Ref. starts with FeAs, Gd, Gd 2 3 and GdF 3 in 
a stoichiometric ratio. All materials were homogenized 
by grinding in a mortar. Route 2 uses GdAs, Fe, Fe2 03 
and FeF3 as starting materials in a stoichiometric ratio. 
Here, the starting materials were homogenized by grind- 
ing in a ball mill. In either case, the resulting powders 
were pressed into pellets under Ar atmosphere, and sub- 
sequently annealed in an evacuated quartz tube either in 
a two step synthesis at 940°C for 12 h and at 1150°C 
for 48 h (60 h) or in a one step synthesis at 940° C for 
168 h. In order to confirm the single phase character of 
the polycrystals, powder x-ray diffraction was performed 
on a Rigaku diffractometer (Cu A Q -radiation, graphite 
monochromator) . The samples were either phase pure or 
contained insignificantly small amounts of GdAs, GdOF, 
and Fe3 04. The microstructure and the composition 
were examined by scanning electron microscopy (SEM, 
XL30 Philipps, IN400) equipped with an electron mi- 
croprobc analyzer for semi-quantitative elemental analy- 
sis using the wave length dispersive x-ray (WDX) mode. 
The analysis showed that the sample with x = 0.15 (nom- 
inal content) in fact contains ~ 0.07 ± 0.02 of F and 
x = 0.17 contains ~ 0.14 ± 0.02 of F. Further on, we will 
use the doping levels obtained by WDX in order to label 
the samples. 

For the c-axis alignment of the parent GdOFeAs sam- 
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FIG. 1: Powder x-ray diffraction data of the c-axis aligned 
GdOFeAs sample. 



pic, which was synthesized under high pressure^, the 
powder was mixed with epoxy resin and hardened while 
rotating in a magnetic field of 1.5 T. The x-ray diffrac- 
tion data of the aligned powder samples were collected 
at room temperature using a PANalytical X'Pert PRO 
system (Philips) with Co A a -radiation (Fig. [J). The 
presence of highly intense [001] reflections (Fig.[IJ ar- 
rows) which dominate the pattern points to a sufficiently 
good quality of the alignment. Reflections with Miller 
indices different from [001] (Fig.[TJ asterisks) are visible 
in the background, too, but their intensity is strongly 
suppressed compared to the powder pattern. 



III. THERMODYNAMIC AND TRANSPORT 
MEASUREMENTS 

In Fig.[2ja), the temperature dependence of the static 
susceptibility x = M/B of GdOFeAs is presented. As 
x(T) is dominated by the response of the Gd moments, 
the results are very similar for the F-doped samples, 
which are not shown. In general, the data obey the 
Curie- Weiss law which is expected due to the presence 
of paramagnetic Gd 3+ ions. Note that the response of 
the FeAs-layers which is e.g. visible in LaOFeAs is about 
3 orders of magnitude smaller and hence masked by the 
magnetism of the rare earth ions.— The linear temper- 
ature dependence of the inverse susceptibility demon- 
strates the Curie- Weiss-like behavior. Analyzing the data 
in terms of the Curie- Weiss-law yields the antiferromag- 
netic Weiss temperature 6 = — 16 ± IK and the effective 
magnetic moments p e g = 7.81 ±0.04 /xb which is close to 
the magnetic moment of a free Gd 3+ ion (p c ff = 7.94 ^b)- 
At a low temperature of about ~ 5 K there is a kink of 
the magnetization due to the AFM ordering of the Gd 
moments. 

While the structural and Fe magnetic phase transitions 
are not visible in the magnetization data, there are pro- 
nounced anomalies in the specific heat c p and the thermal 
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FIG. 2: (Color online) (a) Static susceptibility \ = M / B 
(left axis), x™ 1 (right axis); (b) specific heat c p (left axis) 
and thermal expansion coefficient a (right axis) of GdOFeAs 
vs. temperature; specific heat data for x = 0.07 and x = 0.14 
samples is artificially shifted down for clarity; inset in (b) 
panel presents specific heat anomaly associated with long- 
range antiferromagnetic ordering of the Gd-moments at B — 
and B = 9T. 



expansion coefficient a (Fig.^b)) in the case of the par- 
ent GdOFeAs sample. There is one broad feature visible 
in the specific heat data. In contrast, the thermal expan- 
sion coefficient exhibits two huge anomalies with opposite 
sign which can be attributed to the structural and SDW 
transitions of the compound at T$t = 136 ± 5 K and 
Tsdw = 128 ± 2 K. In addition, the specific heat data 
reveal a sharp anomaly at T N Gd = 3.8 K which is as- 
sociated with the onset of long range antiferromagnetic 
order of the Gd moments, in accord with the magnetiza- 
tion data. Note that the anomaly is not present in our 
thermal expansion data due to the restricted tempera- 
ture range T > 6K. Upon application of external mag- 
netic fields, Gd order is strongly suppressed as shown 
in Fig. [2] (inset in the lower panel). T N Gd is shifted to 
2.5 K in an external magnetic field of B = 9T. While 
anomalies associated with Gd-ordering are still observed 
in the specific heat data of the F-doped samples with 
x = 0.07 and x = 0.14, there are no visible anomalies 
at higher temperatures (Fig.^b)). This evidences the 
absence of long-range SDW order in the doped, super- 
conducting samples. 

Fig.[3]^a) shows the temperature dependence of the 
electrical resistivity p of the GdOi_ x F x FeAs samples for 
all three doping levels: x = 0, 0.07 and 0.14. To get 
a better insight into the data we present the tempera- 
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FIG. 3: (Color online) (a) Electrical resistivity p of 
GdOi-^F^FeAs samples, x = 0, 0.07, 0.14; (b) electrical re- 
sistivity derivative dp/dT for x — (left axis), x = 0.07, 0.14 
(right axis). 

ture derivatives dp/dT in the bottom panel (Fig.[3](b)). 
The resistivity of the undoped material exhibits features 
closely connected to the structural and magnetic phase 
transitions: a maximum close to T$t and an inflec- 
tion point at Tsdw = 128 K, which are characteristic 
for all the 1111 parent compounds^ ' 17 ' 20 . With dop- 
ing, the electrical resistivity drastically changes its be- 
havior, superconductivity emerges at low temperature 
and the intermediate temperature maximum disappears. 
The SC temperatures T c for x = 0.07 and x = 0.14 
samples amount to 20 K and 45 K, respectively. No pro- 
nounced features of the SDW phase are present in these 
compounds in the whole investigated temperature range. 
However, the normal state behavior of p(T) for x = 0.07 
is very unusual. At high temperatures the resistivity is 
linear down to approx 200 K, then it develops a cur- 
vature and drops below the linear approximation of the 
high temperature part. With decreasing temperature fur- 
ther, p(T) becomes linear again and develops a slight 
opposite curvature at T < 50 K, prior to the onset of 
superconductivity. Upon increasing the F doping level 
in the samples, namely to x = 0.14, this anomaly be- 
comes weaker. A similar drop of p(T) at T < 200 K 
as found here has previously been observed for other 
1 Ill-type pnictide superconductors^ ' 17 ' 21 as well as for 
Bai_ a K g Fe2As 2 22 ' 23 . The qualitative resemblance to the 
sharp drop at Tsdw which is observed in the respective 
parent compounds suggests that the resistivity drop in 
the superconducting samples is indicative of remnants of 
the SDW phase. In fact, a recent study of the Nernst 
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effect on LaOi-^F^FeAs provides strong evidence that 
precursors of the SDW phase develop in the vicinity 
of the resistivity anomaly despite the absence of static 
magnetism 2 ^. 



IV. ELECTRON SPIN RESONANCE, RESULTS 
AND DISCUSSION 



A. GdOFeAs, 9.6 GHz X-band measurements 

The ESR measurements performed at a frequency of 
9.6 GHz on the c-axis aligned GdOFeAs sample in the 
whole temperature range of study and for both sample 
orientations (H ext | or _L c) reveal one broad line with 
the g- factor ~ 2 (Fig.[4|a), inset). Such ESR response 
is typical for the systems where Gd 3+ ions occupy reg- 
ular positions in the crystal lattice with short distances 
between neighboring ions^ 5 -. The Gd 3+ is an S-state ion 
with a half-filled 4/ shell which yields an isotropic g- 
factor equal to 2 and a spin value of 7/2. The rather 
big spin value leads to strong magnetic dipolc-dipole in- 
teractions which together with the unresolved fine struc- 
ture broaden the ESR line. This broadening mechanism 
should lead to a gaussian line shape which is however 
not observed in the spectra. Instead, the lorentzian 3 ^ 
function had to be used to fit the spectra (thin line in 
the inset in Fig. Ufa)) in order to obtain accurate val- 
ues of the resonance field (H res ) and the linewidth. The 
lorentzian shape suggests that homogeneous narrowing of 
the line does take place, which can be caused by isotropic 
exchange interaction between Gd spin a 25 ' 26 . The temper- 
ature dependencies of H res and the linewidth are shown 
on Fig.@Ia),(b). With lowering the temperature no dras- 
tic changes are seen in H res down to ~ 10 K where there 
is a strong shift of the line due to the ordering of Gd 
moments (Fig.@Ja)). The linewidth, in contrast to the 
resonance field, shows clear change in the behavior at 
Tsdw = 128 K for both sample orientations (Fig.[4jb)). 
At temperatures above Ts dw there is a gradual decrease 
of the width of the ESR line upon cooling. This can be 
attributed to a Korringa-like behavior, with the linewidth 
having a linear in T contribution due to the relaxation 
of the Gd spins through interaction with the conduc- 
tion electrons, similar to EuFe2As2^1. The slope value 
amounts to ~ 0.9 • 10~ 4 T/K which is one order of mag- 
nitude smaller than that in EUFC2AS2. A strong broad- 
ening of the line below Tsdw can be attributed to the 
formation of the SDW state in the FeAs layers. Similar 
effects in the Gd 3+ ESR linewidth were observed before 
in the case of Gd2BaCu0.5 samples where exchange cou- 
pling of Gd- and Cu-moments enabled to probe by means 
of Gd 3+ ESR the magnetic ordering of the Cu layers^ 
(see the discussion below). 




■g 0.0' 0.2 0.4 0.6 0.8 1.0 




□ H||c 



i , - - ~ v y 4'--'-, -," T -^w w y — _ n , — — t 

50 100 150 200 250 300 

Temperature (K) 

FIG. 4: (Color online) Results of the X-Band ESR (u = 
9.6 GHz) measurements performed on c-axis aligned powder 
GdOFeAs sample for two sample orientations in the magnetic 
field, open circles - field in ab plane, open squares - field par- 
allel to the c-axis; (a) temperature dependence of the res- 
onance field on a reduced field scale (H — Ho) /Ho- Here 
fMoHo = 0.323/0.328 T is resonance field at the highest tem- 
perature for H _L c/H || c; the inset shows the spectrum at 
T — 70 K in H || c configuration, the arrow points at a small 
Fe ESR signal which presence indicates a small amount of im- 
purities; (b) temperature dependence of the linewidth; inset 
shows the change of the T-dependence at the SDW transition 



B. GdOFeAs, high- frequency /field measurements 

In the measurements performed at 9.6 GHz, the 
linewidth of the Gd 3+ ESR signal is comparable to 
its resonance field. This leads to complications in the 
spectra analysis and to a lack of resolution. In or- 
der to improve the spectral resolution we performed 
high-frequency/field measurements on GdOFeAs sam- 
ples. The high-temperature ESR spectra of the non- 
oriented GdOFeAs powder sample measured at a fre- 
quency of 328 GHz (Fig.[5[a)) consist of a single broad 
lorentzian-shaped line with a g-factor of ~ 2 and a 
linewidth of ~ 0.2 T, similar to the low-frequency mea- 
surements. However, the very small Korringa con- 
tribution detected in the low frequency measurements 
is not visible in the HF-ESR spectra (Fig.EK). The 
low-temperature HF-ESR spectra exhibit an inhomoge- 
neously broadened shape which is in contrast to X-band 
data. As a measure of this broadening the full width at 



5 




(H-H )/H (%) (H-H )/H (%) 

FIG. 5: (Color online) Temperature evolution of the high- 
frequency /field ESR spectra of GdOFeAs powder and c-axis 
aligned powder samples at a frequency of v — 328 GHz, shown 
on a reduced field scale (H — Ho) /Ho- Here fXoHo = 11.7 T 
is the resonance field of the signal at high temperature; (a) 
non-oriented powder; (b) c-axis oriented powder. 
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FIG. 6: (Color online) The shift of the minimum of absorption 
of the spectra with temperature, measured at a frequency of 
v = 328 GHz, on a reduced field scale (H rea — Ho)/ Ho- Here 
Ho Ho = 11.7 T is the resonance field of the signal at high 
temperature; open circles - powder sample, open squares - c- 
axis aligned powder sample, solid lines - fits of the resonance 
field using Eq. IA.4I 



the half maximum (FWHM) of the signal AH has been 
taken (see Fig.|SJa)). As can be seen, with decreasing 
the temperature there is only a weak broadening of the 
signal down to a characteristic temperature Tif ~ 150 K 
(if denotes an internal field at the Gd ion, see the discus- 
sion below) where substantial inhomogencous broadening 
begins to develop continuously down to the lowest mea- 
sured temperature of 4K. Concomitantly with the inho- 
mogeneous broadening there is a noticeable shift of the 
minimum of the absorption to lower fields, as shown on 
Fig.[6](open circles) on a reduced field scale (H—Hq)/Hq. 
Here hqHq = 11.7 T is the resonance field of the signal 
at high temperature. The spectral shape at low tem- 



perature (Fig.[5Ja)) appears to be very similar to the 
shape of the ESR signal from a powder sample with an 
anisotropic (/-factor. However, since the Gd 3+ ion is a 
pure S-state ion, it should have an isotropic (/-factor very 
close to 2. Hence one can conjecture that the shape of the 
ESR signal from the GdOFeAs sample is caused by the 
anisotropy of the internal field at the Gd site arising from 
the AFM-ordered Fe moments. In such an anisotropic 
powder situation most of the spectral weight is coming 
from the grains whose c-axes are oriented perpendicular 
to the direction of the external field H cxt . One should 
note here that in the case of in-plane anisotropy there will 
be an additional averaging effect due to the distribution 
of resonance fields of grains whose (a6)-planes are parallel 
to -f/ext ■ Therefore we assume here that the low-field min- 
imum of the absorption corresponds to the mean value 
of the resonance field of the Gd 3+ ESR response (H rea ±) 
in the case of the external field applied perpendicular to 
the c-axis. Correspondingly, the high-field shoulder of 
the spectra arises from grains whose c-axes make small 
angles with respect to H eXt . 

In order to probe the Gd 3+ response for the geometry 
-Hext || c we have performed ESR measurements on the c- 
axis oriented GdOFeAs powder sample. Though, accord- 
ing to the x-ray diffraction analysis, the alignment of the 
powder particles was not perfect, a substantial c-axis tex- 
turing of the sample has been achieved (Fig.[T]). Similarly 
to the non-oriented powder sample, the c-axis oriented 
sample at temperatures above ~ 150 K exhibits only a 
small broadening of the ESR spectrum with decreasing 
temperature (Fig.^b)). Below this temperature the sig- 
nal experiences strong inhomogencous broadening where 
most of the spectral weight is shifted to higher fields, 
which is opposite to the finding in the non-oriented pow- 
der sample. In the oriented sample most of the spectral 
weight and, consequently, the minimum of the absorp- 
tion should correspond to the resonance field of the Gd 3+ 
ESR response (H res \\) in the geometry H ext \\ c whereas a 
non-ideal powder alignment yields the low-field shoulder 
of the ESR signal. 

From our measurements on non-oriented and c-axis 
oriented powder samples one can, therefore, extract the 
temperature dependencies of H res in two configurations, 
i.e., for fields aligned along the c-axis (H res \\) and in 
the (afr)-plane (H res ±), as summarized in Fig. [5] As can 
be seen, the changes of both resonance fields H res ± and 
H res ii start upon cooling at T%f < 150 K and the shifts 
have opposite directions. 

The qualitative difference in the high- field/frequency 
and low-field/frcqucncy measurements leads to the con- 
clusion that the shift of the resonance field and the in- 
homogencous broadening of the spectra measured at a 
frequency of 328 GHz is a field-induced effect. To in- 
vestigate it, we have measured the frequency v versus 
magnetic field H les dependence of the GdOFeAs pow- 
der and the c-axis aligned samples, respectively, both at 
T = 280 K and T = 4K (Fig-.Q. At T = 280 K the 
spectrum at all studied frequencies and fields consists of 
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FIG. 7: (Color online) Frequency dependence of the ESR sig- 
nal from GdOFeAs non-oriented powder and c-axis oriented 
powder samples measured at T — 4K; data points correspond 
to the position of the minimum of absorption of the spectra at 
different frequencies, open circles - non-oriented powder sam- 
ple, open squares - c-axis oriented powder sample, dashed 
lines - guides for the eyes, solid line - the frequency depen- 
dence of the position of the high-T Gd ESR line measured 
at T = 280 K; the lineshape of the spectra is shown as well, 
thin solid line - non-oriented powder, thick solid line - c-axis 
oriented powder. 

a single lorentzian line with the same linewidth value. A 
linear v(H res ) dependence has been revealed yielding a 
^-factor g = hv j (fj,BH res ) equal to 2.005 (straight solid 
line on Fig. [7]). The spectra at T = 4K for the non- 
oriented powder (thin line) and c-axis oriented powder 
(thick line) samples at different frequencies together with 
the frequency dependence of the resonance fields H res ± 
and H res u are shown on Fig.[7]as well. This measurement 
reveals that the difference between H res ± and H res u in- 
creases linearly with increasing the frequency and the 
field strength. 

C. GdOFeAs, discussion 

Since the inhomogeneous broadening and the shift of 
the Gd ESR signal set in close to the SDW ordering tem- 
perature, one can associate them with the interaction of 
the Fe ordered moments with the Gd spins. The shifts 
of H res ± and H res \\ from a common high-temperature 
paramagnetic value Hq can hence be related to the oc- 
currence of an internal magnetic field (Hi nt j_ and H int u) 
at the Gd site due to the formation of the static SDW in 
the FeAs layer. Based on the dependence shown on Fig. [7] 
we can conclude that the strength of these internal fields 
depends on the strength of the applied magnetic field. 

When the FeAs planes are in the paramagnetic 
state, at temperatures above Tsdw = 128 K, the ex- 
change/dipolar field at the Gd position is negligible as the 
applied magnetic field cannot effectively polarize small 



paramagnetic moments at elevated temperatures. At 
temperatures below Tsdw = 128 K the Fe moments or- 
der statically in the (a&)-plane. The AFM structure of 
the magnetic order, which is similar for all i?-based 1111- 
type pnictides can be found in Ref. HH Due to symmetry 
reasons, without an applied magnetic field the internal 
field at the Gd site can be only of dipolar nature and in 
this case according to an estimate it should not exceed 
~ 0.03 T. We suggest that application of an external mag- 
netic field induces tilting of the Fe spins and hence creates 
an uncompensated magnetic moment (mp e ) in the direc- 
tion of the external field^ 9 -. This moment can interact 
with the Gd 3+ spins which would effectively lead to the 
occurrence of an additional internal field. If H ext _L c, 
then the Fe spins tilt in the ab plane as it is schemat- 
ically shown on Fig.(5][a,b) for the configuration where 
H ext makes an angle of 90° with the AFM ordered Fe 
spins. The shift of the ESR line H res ±, measured at 
328 GHz, yields an estimate of the internal field (Hi nt ±) 
of about ps 0.4 T parallel both to H ext and to the un- 
compensated moment mp e in the FeAs plane. If H ext \\ c, 
then the Fe moments tilt out of the plane (Fig. Etc, d)). In 
this case the shift of the ESR line, measured at 328 GHz, 
yields an estimate of the internal field (H int u) of about 
w 0.65 T antiparallel both to the H ext and to the uncom- 
pensated moment mp e . One should note here that an 
estimate of the dipolar field produced by the FeAs lay- 
ers at the Gd site yields a value which does not exceed 
pa 0.05 T even for full out-of-plane canting. This field is 
one order of magnitude smaller than the experimentally 
observed value which clearly implies the presence of an 
appreciable exchange interaction between Gd and canted 
Fe moments. The dependence of the sign of the internal 
field on the direction of the applied magnetic field sug- 
gests that the sign of the exchange interaction with Gd 
spins is different for different directions of the uncompen- 
sated Fe moments, i.e., ferromagnetic for the in-plane and 
antiferromagnetic for the out-of-plane directions. This 
surprisingly strong anisotropy of the exchange might be 
related with the multiband electronic structure of iron 
pnictides which might give rise to different pathways for 
interactions between the Gd 4/ orbitals and the in-plane 
xy and out-of-plane xz and yz Fe 3d bands. Note that, 
in zero magnetic field and hence without Fe spin canting 
(rrip e = 0), the exchange interaction between Fe and Gd 
moments of an arbitrary sign is geometrically frustrated 
(see Fig. [8]). The application of a field which tilts the Fe 
moments thus removes this frustration. 

Considering the exchange interaction, one might sup- 
pose that the temperature dependence of the internal 
field at the Gd site should follow the behavior of the 
SDW order parameter in 1111 compounds ^ 15 ' 29 , which 
increases fast within ~ 30 K starting at Tsdw and then 
stays almost constant with further decreasing the tem- 
perature. Here, the internal field acting on the Gd 3+ 
moments arises at a temperature Tif ~ 150 K which is 
~ 20 K higher than Tsdw = 128 K and keeps increas- 
ing upon cooling till the lowest measured temperature 
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FIG. 8: (Color online) Canting of the Fe moments due to the 
applied magnetic field H ex t', (a),(b) H ex tA- c (angle between 
H ext and Fe AFM ordered moments is 90°); (c),(d) H ex t \\ 
c; arrows on the Fe site depict magnetic moments, whereas 
arrows on the Gd sites represent the induced internal field. 



(see Fig.El) • When approaching the temperature of the 
SDW transition from high temperatures, the appearance 
of the internal field well above Tg dw can be explained 
by growing quasi-static correlations between the Fe mo- 
ments seen in the time window of the high-frequency 
ESR. The development of the internal field below T$dw 
is found to be similar to that of some other systems where 
paramagnetic ions are coupled to magnetically-ordered 
moments of another typ o 30 ' 31 . To explain this evolution 
of the internal field we use a simple model based on a 
mean-field approximatio n 25 ' 32 (for details see Appendix). 
According to this model the internal field depends on the 
applied magnetic field and on temperature as following 



Hi, 



± 



A(J,a,C Gd )H e 
(T-6) 



Here A{ J, a, Cad) is a parameter determined by the Gd- 
Fe exchange coupling energy J, by the susceptibility of 
the ordered Fe moments to the applied field a and by 
the Gd Curie constant Cad- As can be seen, this depen- 
dence qualitatively obeys a Curie- Weiss law at a given 
applied field which agrees well with the measured data. 
The model enables to fit the experimental data points 
reasonably well (see Appendix and Fig. [6]). Firstly, the 
fit yields an estimate of the energy of the exchange cou- 
pling between Gd and Fe spins which is in the range 
of |J| w 15 — 20 K. Secondly, it shows that the magni- 
tude of the uncompensated Fe moment depends on the 
Gd magnetization which suggests that the Gd subsystem 
additionally tilts or polarizes the SDW. 
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FIG. 9: (Color online) The full width at the half maximum 
(FWHM) AH of the ESR lines, measured at a frequency of 
v = 348 GHz (x = 0.07) and v = 328 GHz (x = 0,0.14), as a 
function of temperature for the non-oriented GdOi-iF^FeAs 
samples; (a) x = 0; (b) x = 0.07; (c) x = 0.14. 



To summarize this part, the ESR results on GdOFeAs 
samples show that the Gd subsystem is exchange-coupled 
to the magnetic FeAs planes. On approaching the AFM 
SDW transition from above, the growing correlations be- 
tween the Fe moments yield a shift of the Gd ESR line. 
At lower temperature, depending on the angle between 
H ext and the c-axis of the sample, the signal shifts to 
higher or to lower fields due to the uncompensated ex- 
change field which is transferred to the Gd site from the 
Fe moments canted in an external magnetic field. Since 
the full width at the half maximum AH is proportional 
to the difference between resonance fields H res ± — H res \\ , 
then the width of the ESR signal AH of the non-oriented 
powder sample can be taken as a measure of this ex- 
change field (Fig.[9K). 



D. 



GdOi-^FeAs (x = 0.07, 0.14) 



The influence of the fluorine doping on the Gd ESR has 
been studied on two powder samples of GdOi-^F^FeAs 
with 7% and 14% of fluorine. On Fig-fTOfcUb) the evo- 
lution of the respective Gd 3+ ESR spectra is shown on 
a reduced field scale. Similarly to the undoped sample, 
at high temperature the ESR spectrum for both doped 
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(H-H o )/H (%) (H-H )/H (%) 

FIG. 10: (Color online) Temperature evolution of the high- 
frequency /field ESR spectra of GdOi-^F^FeAs powder sam- 
ples measured at a frequency of v — 348 GHz (x — 0.07) 
and v = 328 GHz (x — 0.14), shown on a reduced field scale 
(H — Hq) /Ho. Here Ho is the resonance field of the signal at 
high temperature; (a) x = 0.07, n H = 12.4 T; (b) x = 0.14, 
lioHo = 11.7 T. 




Temperature (K) 



FIG. 11: (Color online) The shift of the minimum of absorp- 
tion of the spectra with temperature, measured at a frequency 
v = 348 GHz(a; = 0.07) /328 GHz (a; = 0.14), on a reduced field 
scale (H res - H )/H Q . Here fi H = 12.4T/11.7T is the res- 
onance field of the signal at high temperature, open hexagons 
- x = 0.07, open triangles - x = 0.14, solid line - fit of the 
resonance field in the parent GdOFeAs sample. 



samples consists of a single lorentzian-shaped line with 
g = 2.005. While lowering the temperature, the line re- 
mains almost unchanged until a characteristic tempera- 
ture Ti f is reached. This temperature corresponds to the 
onset of an additional inhomogeneous contribution to the 
width of the ESR signal AH. This inhomogeneous con- 
tribution is shown by the shaded area on Fig.[S£b,c). The 
temperature Tjj clearly depends on the fluorine doping 
level. In the case of the 7% F-doped sample a noticeable 
broadening of the line starts at T{f(x = 0.07) ~ 125 K, 
whereas for the 14% doped sample it starts at a lower 
temperature Ti/(x = 0.14) ~ 100K. For both doped 




FIG. 12: (Color online) Frequency dependence of the ESR 
signal from GdOi-aT^FeAs powder samples measured at a 
temperature of T = 4K, in SC state; data points correspond 
to the position of the minimum of absorption of the spectra at 
different frequencies, open hexagons - x = 0.07, open triangles 
- x = 0.14, dashed line - guide for the eyes, solid line - the 
frequency dependence of the position of the high-T Gd ESR 
line measured at T = 280 K for both F-doped samples; the 
lineshape of the spectra is shown as well, thin solid line - 
x = 0.07, thick solid line - x = 0.14. 



samples there is a shift of the minimum of the absorp- 
tion (H res ±) to lower magnetic fields below Tif (Fig.lTTl). 
Qualitatively this shift is similar to that of the undoped 
sample (solid line on Fig. lll|) . but it is less pronounced. 
In addition, the inhomogeneous broadening and the shift 
of H res ± to lower magnetic fields exhibit a magnetic field 
dependence similar to that of the undoped GdOFeAs 
sample (Fig. [T2"|) . 



E. Superconducting GdOi ^F^FeAs, discussion 

The remarkable similarities of the Gd ESR behavior 
between the fluorine-doped samples and the undoped one 
strongly suggest that, even in the superconducting sam- 
ples where the phase transition to the AFM SDW state is 
not observed in the thermodynamics and transport prop- 
erties, quasi static (on the time scale of the ESR mea- 
surement) magnetic correlations in the FeAs planes are 
present below the characteristic temperature Tif. Such 
correlations may explain the peculiar features in the re- 
sistivity data shown in Sec. IIIII 

The unified phase diagram for the iron pnictides, es- 
pecially for the 1111 materials, is not fully established 
so far since the issue of coexistence of superconductivity 
and magnetism remains controversial. The nonmagnetic 
rare earth based LaOi-^F^FeAs material exhibits no ev- 
idence for the presence of a static magnetic order for any 
superconducting composition, but rather reveals SDW- 
likc spin fluctuations seen in the transport^! and inelastic 
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neutron scattering experiment o ' . On the other hand, 
the magnetic rare earth based systems (Sm, Nd, Ce) 
studied so far demonstrate coexistence of superconduc- 
tivity and static magnetism at least in the underdoped 



region 3 -^. Our high-field ESR results show that there 
is yet another 1111 system comprising a strongly mag- 
netic rare earth (Gd) subsystem where the coexistence 
of quasi-static magnetism and superconductivity is still 
visible in large doping range. Here, the increase of the 
fluorine content and correspondingly the rise of T c leads 
to the suppression of magnetic correlations indicating a 
possible interplay between these two states. All this sug- 
gests that the coexistence and possible interplay of the 
static or quasi-static magnetism and superconductivity 
may be a generic property of 1111-type compounds. In 
this regard, a remaining question yet to be answered is 
the extent to which the R-Fe magnetic interaction in- 
fluences the magnetic correlations in the FeAs planes. 
Moreover, the issue of the coexistence of magnetism and 
superconductivity is frequently discussed in the literature 
for other pnictidc families as well. Extending ESR exper- 
iments to compounds of these families would be of great 
interest. 



V. CONCLUSION 

Our HF-ESR study of polycrystalline samples of the 
GdOi_ x F x FeAs superconductor reveals a magnetic cou- 
pling between the Gd subsystem and the FeAs layers. 
This coupling, most probably of the anisotropic exchange 
type, is visible in the Gd ESR response in the undoped 
GdOFcAs in the SDW state, in form of a field-induced in- 
homogeneous broadening and shift of the ESR spectrum. 
This effect is caused by the interaction of the Gd spins 
with the uncompensated Fe moments due to the canting 
of the Fe moments in magnetic field. Furthermore, the 
data suggest that the Gd moments additionally tilt the 
ordered Fe moments. Surprisingly, the broadening and 
the shift of the spectrum are present also in the doped 
superconducting samples where there is no evidence of 
long range magnetic order. This points to the presence of 
short range, static on the ESR time scale, magnetic cor- 
relations. This may be relevant to the interplay of mag- 
netism and superconductivity in these materials, where 
on doping with fluorine there is a simultaneous increase 
of the superconducting critical temperature and suppres- 
sion of the magnetic correlations. The possible relevance 
of the exchange interaction between the magnetic rare- 
earth subsystem and the FeAs planes to the properties 
of this novel class of superconductors remains to be elu- 
cidated. 
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Appendix 

Here we provide a calculation of the internal field on 
the Gd ion using a simple model based on the mean-field 
approximatio n 25 ' 32 . To simplify the calculations we as- 
sume that a magnetic field and magnetic moment vectors 
are collinear. In this model, the internal field at the Gd 
site Hi n t is proportional to the magnitude of the uncom- 
pensated Fe moment mp e with a coefficient A: 



Hi, 



±Xm Fe 



(A.l) 



Hereafter the sign depends on the type of interaction, be- 
ing "+" for ferromagnetic and "— " for antifcrromagnetic 
exchange. Neglecting the weak dipolar contribution, the 
magnetization normalized to the single ion mad of the 
Gd subsystem is proportional to the sum of the applied 
magnetic field H ext and internal field Hint- 



m G d = XGd(T) (H ext ± Hint), 



(A.2) 



where XGd{T) = Cad/iT — 0) is the Gd magnetic sus- 
ceptibility, Cad is the Gd Curie constant, O is the Gd 
Curie temperature. Due to the exchange interaction, the 
uncompensated magnetic moment m Fe is proportional 
not only to the applied field but also to the internal field 
created by the Gd moments: 



m Fe = a{H ext + Xm G d) 



(A.3) 



Here a is the susceptibility of the ordered Fe moments to 
the external magnetic field. Using Eg. IA.11 Eg. IA.2l and 
Eg. lA.3l one can obtain the equation for the internal field 
Hi n t at the Gd site: 



Hint — i 



X 2 aC Gd + Xa{T - O) 
(T — 9) - X 2 aC Gd 



H 



ex I 



(A.4) 



Eq. IA.4I enables to fit the measured temperature de- 
pendence of the internal field (see Fig-EJ) - The resonance 
field of the Gd is determined by the applied field and by 
the internal field Hq = H ext ± Hint- At high tempera- 
tures when there is no internal field at the Gd site at any 
strength of the applied magnetic field {Hint = 0) one can 
measure the resonance field Hq = 11.7 T (for measure- 
ment frequency v = 328 GHz). Assuming that the reso- 
nance field of the Gd ions Hq = 11.7 T stays constant at 
all measured temperatures one obtains the expression for 
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the applied field H ext = H ^H int (H = U.7T). The fit 
for two measurement configurations {H res ± and -ff re s||) 
is shown on Fig.|6]by solid lines. The parameters Cad-, @ 
and a can be taken from different experiments. The Gd 
Curie constant C'cd and Curie temperature are known 
from the susceptibility data of GdOFeAs samples (see 
Sec. IIII[) . As it is shown in Ref. [l!^ the bulk Fe suscep- 
tibility of LaOFcAs samples is determined by the spin 
susceptibility. Therefore the parameter a can be esti- 
mated from this measurement yielding a value of ~ 10~ 4 
. The A value resulting from the fit is equal 
to ~ 19.7 ^Jfi for H res± and ~ 25.3 for H rest 

According to the mean field theory^, these values yield 



an estimate of the exchange interaction energy J for two 
configurations amounting to |J| ~ 15 K for H res ± and 
|J| ~ 19 K for H res \\. In addition, Eq. [AT] enables to cal- 
culate the uncompensated moment mpe- Its value grows 
with decreasing the temperature and increasing the Gd 
susceptibility until it reaches ~ 0.03/is at the lowest mea- 
sured temperature. 

Since a is very small compared to Cad and A, Eq. IA.4I 
can be simplified to the form: 



Hi, 



^ X 2 aC Gd 
± (T-Q) Mext 



(A.5) 
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